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Abstract 

An exact calculation of the eigenvalue statistics of truncated random 
Haar distributed real orthogonal matrices has recently been carried 
out by Khoruzhenko, Sommers and Zyczkowski. We further develop 
this calculation, and use it to deduce a Pfafnan form of the correlations 
for the zeros of the limiting Kac random polynomial. This contrasts 
with the forms known from previous studies of the real zeros (a mul- 
tidimensional Gaussian integral with the integrand multiplied by the 
absolute values of the variables) and the complex zeros (a Hafnian). 



1 Introduction 

There has been much recent progress in the exact evaluation of correlation func- 
tions for various non-Hermitian, real entried random matrix ensembles [U [9j \19\ 
EQl \7\ El [U [13] . The most recent of these works, due to Khoruzhenko et al. [13], lays 
the foundations to completely solve this problem for the ensemble of real entried 
random matrices formed from M x M sub-blocks of (L + M) x (L + M) Haar 
distributed real orthogonal matrices. In light of the recent work of Krishnapur 
|14j . the case L = 1, M — > oo, has immediate implication to the evaluation of the 
N — > oo limiting correlation functions for the zeros of the so-called Kac random 
polynomial 

N 

f{x) = Y,<h*?> (!) 

p=0 

where each a p is an independent standard real Gaussian. It is the purpose of 
this paper to make these implications explicit, and to similarly make explicit the 
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implications for the point process denned by the zeros of the matrix-valued real 
Gaussian power series 

oo 

det(^^ p) x p ) =0. (2) 

p=0 

Here each is an independent L x L standard Gaussian real random matrix. 

We begin in Section 2 by adapting the work of Krishnapur [14] to show the 
relationship between the eigenvalues of truncations of random orthogonal matrices, 
and the zeros of ([2]). We further revise known exact expressions, due to Bleher and 
Di [3], for the correlations of the real zeros of (pQ), and known exact expressions 
due to Prosen [18] for the correlations of the complex zeros. In Section 3 we 
extend the workings in [13] in relation to the exact evaluation of the correlations 
for eigenvalues of truncations of random orthogonal matrices. Thus we give the 
explicit M — > oo form of the entries of the Pfafnan specifying this correlation 
function. 

The result from Section 2 giving an equivalence between the eigenvalues of 
certain truncated orthogonal matrices and the limiting zeros of ([2D, together with 
the Pfafhan form of the eigenvalue correlations obtained in Section 3 tell us that 
the limiting zeros of the Kac polynomial form a Pfafnan point process. However, 
neither the correlations between the real zeros, nor the correlations between the 
complex zeros obtained in [3] and [18] respectively, are in a Pfafhan form. We 
can check directly that the one and two point function of the former, and the one 
point function of the latter, coincide. This is the essence of our conclusion, Section 
4, which therefore leaves open the problem of directly showing the equivalence 
between the different forms. 



2 Preliminary results 
2.1 The result of Krishnapur 

Consider the ensemble {U} of Haar distributed (L + M) x (L + M) real orthogonal 
matrices. Block decompose U according to 



U 



A LxL B 
C Vmxm 



(3) 



As noted in [14] in the case of complex unitary matrices, these blocks satisfy the 
determinant identity 

det(zI M - V T ) _ / i \M J„ + TT A„i-f A , -DAT.. 



det(I M - zV) 
Now write 



(-1) M det U det(A + zB(I M - zV)~ L C). (4) 



t f\ i i\M A , TT det(zI M -V) 

f M {z) = {-l) detC/— — . 5) 

det(llM — zv) 
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According to (jU 




But we know from [10^ pg. 113] that 



VM(A, BC, BVC, BV 2 C,...) ^(4 0) , Af , Af , . . . ) 



(7) 



and furthermore that for M — > oo this implies the eigenvalues of the sub-block V 
coincide with the zeros of ([2]) for N — > oo (compare ([5]) with ([6]), 0). 

2.2 Limiting real and complex correlations for the Kac 
polynomial 

In relation to ([I]), it was proved by Kac [T2] that for large N the leading order ex- 
pected number of real zeros is (2/ir) log N. Also, it is known (see [3] and references 
therein) that for N — > oo the density of real eigenvalues is 



This is singular at x = ±1 — the points where the zeros accumulate, but is 
otherwise well defined. 

An important point from [5] is that the zeros for |x| < 1 are statistically 
independent from the zeros for |x| > 1. This is in keeping with the fact [10, 
Lemma 2.3.3] that for N — > oo the random analytic function ([2]) has radius of 
convergence unity, and moreover cannot be analytically continued to \x\ > 1.. It 
also explains why it is consistent that the eigenvalues of V in Section 2.1 all have 
modulus less than 1. 

The general /c-point correlation function p T ^ for the real zeros in (—1,1) of the 
limiting Kac polynomial ([T]) has been computed by Bleher and Di [3]. Furthermore, 
these authors have shown that this same functional form persists for a wide class 
of mean zero, variance one distributions for the coefficients a p in JT]) [3]. To state 
the result, let 



Afc — [Aji]j,i=i,...,k 



with each Aji a 2 x 2 block 



1 Xj 




1 - XjXl (1 - XjXl) 2 
Xi 1 + XjXl 



(1 - XjXl) 2 (1 - XjXl) 3 
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Then we have from [H eq. (4.6)] that 
p\ k) {xi,...,x k ) 

roo 

dyi--- I dy k \ yi ---y k \e-y 1Xlk \ (9) 



oo 



(2vr)Vdet A k 

where y = [yi]i=i ,„ } k is a column vector and tt k denotes the k x k matrix obtained 
by removing all the odd number rows and columns from A7 . In the case k = 1 
this reclaims ([8|), while for k = 2 it simplifies to give [U eq. (3.10)] 

r ( \ (Xi -X 2 ) 2 

P{2) (X 1 ,X 2 ) - 7T 2 {1 _ XlX2) 2 {l _ x 2 ){1 _ x 2 ) 



. \xi-x 2 \ . s/{l-x{){\ -sjQm N 

H , arcsm — ^flO) 

7r 2 (l-xia; 2 )V(l-^)( 1 - ;E 2) l-x x x 2 
(for an extended discussion of issues relating to a direct evaluation of the integral 
in Q for general £l k , see [To]). 

We now turn our attention to the complex zeros. For the a p in ([1]) any inde- 
pendent real Gaussian distribution with variance b p , and thus 

(a p a q ) = b p 6 Pi q, (11) 

a formula for the /c-point correlation p c ^ of the complex zeros has been given by 
Prosen [18]. Since the coefficients of the polynomials are real, the complex zeros 
must occur in complex conjugate pairs. The exact evaluation is given in terms of 
a quantity known as a Hafnian (see e.g. [UJ) but referred to as a semi-permanent 
in pg, 

k 

HfA 2 £;x2fc= * JJop(2«-l),P(2Z)- (12) 

Pes 2k 1=1 

Here A 2kx2k = [<ijk]j,k=l ... 2k is assumed symmetric, and the asterisk on the sum 
means that the permutations must statisfy -P(l) < -P(3) < • • • < P(2k — 1) and 
P(2i - 1) < P(2i) (i = 1, . . . , k). According to [Tgl eq. (12)] we have that 

r , ^ Hf(C - B^A- X B) , s 

*> (zi 7mm ■ <13> 

The matrices A,B,C are of size 2k x 2k and have entries defined in terms of 

N 

g(u):=Y,b n u n , (14) 

n=0 

where b n is specified by {UJ). Explicitly, with z k+s := z s and ^'(u) denoting the 
derivative of g(u), these matrices have entries 

Aji = g{zjzi) 

Bji = d Sl Aji = Zjg'(zjZi) 

Cji = d Zj d Zl Aji = g'(zjZi) + zjzig' \zjz~i). 
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In the present problem we have b n = 1 for each n, and we must also take the limit 
N — > oo. Then (j!4j) simplifies to read 



1 



u 



(15) 



In the case k = 1 (fTB"|) can be written out explicitly ([18, eq. (19)] for general 
5. Substituting (fT5|) in that formula we find 



2: — z\ 



vr|l-z 2 |(l 



1 2\2 ' 



(16) 



3 Pfaffian correlations for eigenvalues of trun- 
cated orthogonal matrices 

3.1 Specification of the correlations 

Let an (L + M) x (L + M) random Haar distributed real orthogonal matrix U be 
decomposed as in Q. The exact form of the joint eigenvalue PDF and the /c-point 
correlations has been the subject of the recent work [13] . It was shown that up to 
a known normalization, the joint eigenvalue PDF is given by 



n & 

l<j<fc<M 



M 

Zk)J\w{Zj) 

.7 = 1 



(17) 



where 



w(z) 



f L ( Z, ~ 1 ) | 1 -21L-2 
V 2vr 1 1 

f-) 1/2 ii-^ 2 r 1/2 , 

1 V27r; 1 1 



1 

%i/i 



(l_ u 2 )( L-3)/2 dn \ 
l-z 2 | 7 



1/2 



L = 1. 



(18) 



Note that when 2; is real and thus y = 0, (|T8j) redues to 



w(x) 



l r(i/2)r((L + i)/2)\V 



2tt 



r(L/2) 



1 



,2^/2-1 



(19) 



The functional form ()17p must be considered as a sum over all allowed numbers 
of real eigenvalues. Thus with M even (for convenience) there can be an even num- 
ber of real eigenvalues, while the remaining eigenvalues come in complex conjugate 
pairs. As noted in [13], but not made explicit, the structured form (fT7"|) implies 
that the correlation functions between general sets of k\ real eigenvalues and &2 
complex eigenvalues (the latter with imaginary part positive) can be written as a 
2(fci + k 2 ) x 2{k 1 + k 2 ) Pfaffian 0EIE1 ?]. 
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The key to such formulas is an antisymmetric inner product induced by ([T 



2i 



D+ 
1 r l 



(w(z)f [g 1 {z)g 2 {z) - g 1 (z)g 2 (z))dxdy 



w(xi)w(x 2 )g(xi)g(x 2 )sgn (x 2 - x\) dx\dx 2 (20) 



-1 J -1 



where z = x + iy and D + is the upper half unit disk. 

Let {pj(z)}j=o,i,... he a set of skew orthogonal polynomials with respect to ([20 

(P2j,P2k) = (P2j+l,P2k+l) = 

(P2j,P2k+l) = U^ k ), (j?2j,P2j+l) = Tj, 



(21) 



where it is required that each pj(z) is monic and of degree j. Use these monic 
polynomials and the weight (fTHJ) to define 



w(u) := \2w{u) 
qj(u) := w(u)pj(u) 
1 



Tj(u) 

e(ui,u 2 ) 



2 

iqj (u 

1 

2 
0. 



sgn (« — t)qj(t) dt, u G 



sgn(ui-ii 2 ), Mi,u 2 GM 
otherwise 



(22) 
(23) 

(24) 
(25) 



and 



Af/2-1 



D(x,y) 
I{x,y) 



T i v 



M/2-1 



^ — (q 2j (x)q 2j+1 (y) - q2j+i(a 



(26) 
(27) 



j=0 

M/2-1 



^r( T 2j( x ) T 2j+i(y) - r 2j+ i{x)T 2j (y^j +e(x,y). (28) 

(29) 



3=0 



According to working in [H[T6], the (ki + k 2 )-point correlation function is given 
in terms of a Pfaffian by 



P(fcl+A:2)( x l, ■ ■ ■ iX^; Zi, . . . ,Zk 2 ) 
= Pf 



[K TT {xj,xi)]jj = i,..., kl [K rc (xj,zi)}j=i : ..., kl 

!=1,...,& 2 

[K CJ .{Zj,Xi)]j=r,..., kl [-^cc(^,^)]j=i,-,* 2 

!=l,...,fc 2 !=l,...,fc 2 



(30) 



6 



where 



K CT (z,x) 



Srt(x,y) -D„(x,y) 

I„(x,y) S ir (y,x) 

S CT (z,x) -D cr (z,x) 

I CI (z,x) S IC (x,z) 



K vc (x,z) 



K cc (z 1 ,z 2 ) 



S vc (x,z) -D IC (x,z) 
I TC (x,z) S CI [z,x) 

Sec (zi,z 2 ) -D cc (z 1 ,z 2 ) 
{zi,z 2 ) S cc (z 2 ,zi) 



Here S v , D. ., I v refer to (|26j) - (j39j) : the subscripts indicate if the correspond- 
ing variables are real or complex, which according to (|24p and (|25p effects the 
functional form 

For future reference, we remark that (|30p gives for the one and two point 
correlations of the real eigenvalues 



p\i){x) 
p\ 2) {xi,x 2 ) 



S IV (x, x) 

S TT (x 1 ,x 1 )S TT (x 2 ,x 2 ) 

Srr(xi,X 2 )S„(x 2 ,Xi) + D„ {x\ , X 2 )I Tr (xi , X 2 ] 



(31) 



Also, structurally the same formulas hold for the one and two point correlations 
of the complex eigenvalues — replace the subscripts rr by cc. 



3.2 Skew orthogonal polynomials and summation for- 
mulas 

The practical use of (|30p requires the explicit form of the skew orthogonal polyno- 
mials, as determined by (f20|) and (j2~T|) . To proceed directly in their computation 
appears to be a formidable task. However, as observed in [13] . and in [332 GSH 12] 
for other ensembles of non-Hermitian random matrices, an indirect determination 
is possible. This follows by evaluating the LHS of the identity 

(x - y)(det(x - V {M _ 2)x{M _ 2) )(y - V (M _ 2)x(M _ 2) )) 

M/2-1 

= Yl —{pv+i( x )P2j(y)-Pv+i(y)pv( x ))> ( 32 ) 

3=0 3 

valid for M even. Here Vim—2)x(m— 2) is the bottom right block in ([2D in the case 
that the real orthogonal matrix U is of size {L + M — 2) x (L + M — 2), and the 
average is over Vim-2)x(M-2)- The result of [13], e Q- (8)] is that the LHS is equal 
to 

M-2 , 

Substituting (f33[) in ([32]) . replacing M by M + 2, then subtracting the original 
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equation gives 

— ^— (pM+i{x)pM{y) - pm+i{v)pm{x, 

fM/2 V 

Recalling that each pj(ic) is monic of degree j, and that M is even, it follows from 
this that 

{x^ , j even 

,•' 1 ' ' \ j odd f:5 U 

L + j-1 

and 

I _ (£±M! (35) 
rj " L!(2j)! ' (35) 

Furthermore, comparing the RHS of (|32|) with the definition ([38]) of D(x,y), it 
follows from ([38]) that 

Af-2 

L\m\ 



M~2 , s, 

£>(x, y) = -*(x)u,(y)(x - y) £ 1 -f W". (36) 



m=0 



Recognizing the sum as a partial sum of the generalized binomial series, we obtain 
the summation fomula 



w(x)w(y)(x — y) 
M->oo v ' a ' (1 — xy) 



Hm D( X ,y) = - ^rZT L+ r - (37) 



From knowledge of (|34p we would like to similarly obtain explicit forms for the 
quantities S and / as specified by ([26]) and ([39]) respectively. When both variables 
are complex, this is immediate, as we can check from the definitions that 

Scc(zi,z 2 ) = iD(zi,z 2 ), I cc (zi,z 2 ) = -D(zi, z 2 ). (38) 

Applying ([37]) we thus have 

.w(zi)w{z 2 ){zi - z 2 ) 



lim S cc (zi,z 2 ) - ■ , , , 

M->oc (1 — ZiZ 2 ) L+l 

w(z 1 )w(z 2 )(z 1 - z 2 

Af-»oc LtA X ' Z/ (1 — Z1Z2) 



hm / cc (zi,z 2 ) = ^ - NL+1 . (39) 



Consider next the case that both variables are real. Our strategy to compute 
S TT (x,y) is borrowed from [9]. There structures we are about to exhibit were 
observed in the case of the partially symmetric real Ginibre ensemble. Thus we 
begin by observing from ()34p that 



Recalling the definition (|24j) then shows 



and so 

M/2-1 



Y- 1 , n /n L T(l/2)r((L + 1)/2) 2 ^/ 2 -i 

g = ^ r( L/2) j ( 1 -^ 2 ) L/2 



; (!-y 2 ) L/2 E -T-T^^(x> 2i ( y ). 
j=o 3 J 



Also 



M/2 

E -92j+i(^2j(y) 

i=o ^ 

M/2-1 

J'=0 J J 

= T M -2(y)w(x)x M ~ 1 

TM/2-1 

Mil— 2 

U (L + 2j + l)! 2 . +1 / 2j + l 

g L!(2j + 1)! * (^ +2 ^- L + 2j + l ^ (y) 

To proceed further, we note that analogous to ()40p we have 

«w+^)-r+27+i psw(w) 

= -7 — ! (i - y 2 )" L/2+1 # f (i - y 2 ) L/2 y 2j+1 ) ■ 

L + 2j + V y ' dyV y ' y J 

Recalling now the definition (|24p of Tk we see 

Wv) - TT^TT T ^ y) = L + 2i + i ^ )y2j+1(1 " y2) - 

Substituting (gHJ) in (@2D shows 

M/2-1 x i 
V] — P2j+i(a;)r 2 j(2/) = T M2 (y)w(x)x M 1 

r i r M/2-1 

M/2— 2 

-/ w U1 2n V- (L + 2j + l)! 1 2i+i 

— u;(x)it;(?/)(l — y ) > — — — x J ^ . 

v 1 yy>y y ' p< Q L!(2j + 1)! L + 2j + l 
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Appropriately combining (|4ip and (|44p then gives our sought summation formula 
S tI (x,y) = — T M ^ 2 (y)w{x)x M ~ 1 

r M/2-l 

i r(l/2)r((L + l)/2 _ x2)L/ 2- Ul _ v 2 )L/ 2 A y? + v 
+7T r(L/2) 11 xj ^ (L-i)yi 

Finally, taking the M — > oo limit and relating the resulting sum to the generalized 
binomial series we obtain 

m*. ») - i r<i/2 x; i)/2) ° " T" v y2)L ' 2 ■ <«> 

A/^oo 7T 1 (L/2) (1 — xy) 1 " 

With S rr (x,y) known, I rr (x,y) follows from the formula 

fv i 

I TI (x,y) = - S IT (u,y) du + -sgn (x -y), (46) 



which in turn follows from the definitions (|39|) . (|26|) and (|24|) . We note that in the 
case L = 1 the integral in (|46p with the substitution (|45p can be evaluated to give 



,. 7 f x sgn (y - a) . ^(1 - x 2 )(l - y 2 ) 

nm l ir (x,y) = arcsm . 47 

A/^oo m ' y; l=i vr 1-xy v ; 

Calculation of the quantities S. . and J. . for the real/complex and complex/real 
blocks requires no further work. Thus it follows from the definitions that 

S CT (z,x) = S rr (z,x) 
S TC (x,z) = iD rc (x,z) 

I cx {z,x) = -I IC (x,z) = iS TT (z,x). (48) 



4 Conclusion 



From Section 2.1 we know that the zeros of the matrix- valued real Gaussian power 
series ([2]) are statistically identical to the eigenvalues of the truncated real orthog- 
onal matrices VmxM, M — > oo in (J3J). The correlations for the latter are given by 
the Pfaffian ([21]), with entries specified bv (|37|l . (|55]). ([55 ]) . (jlS]) and ^6] ) . 

In the case L = 1 the matrix valued random function in ([2]) is equal to the 
N — > oo limiting form of the Kac random polynomial ([T]). The explicit form of the 
entries in (13011 are then 



Kn(x,y) 



L=l 



1 

7T 



[i-y 



2U/2 



(x - y) 



(1 -x 2 )V 2 (l - ay) 



sgn(y 



x arcsm ■ 



V(l-X 2 )(l-y2j 

xy 



(l-X 2 ) 1 /2(l_ y 2)l/2( 1 _ xy )2 



X 



2)1/2 



( l_y2)l/2 (1 _ xy) 



(49) 
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K rc (x,z) 



L=l 



-I (X 



(x 



1 

TT 



(1 



jc2)1/2|1. 



2 |V2(1 
,2)1/2 



XZ 



(1 



X 2 )V2|1 



|V2(1 



(1 



_2\l/2 



|l-f 2 |V 2 (l- zx) 



|1/2(1 



(50) 



if cr (^,x) 



L=l 



1 

7T 



(1-x 2 ) 1 / 2 

|1 — J8 2 |V2(1 — zx) 

ijl-x 2 ) 1 / 2 
\l-z 2 \ l / 2 (l-zx) 



— (x — z) 

(1-X 2 ) 1 /2|1_ Z 2|1/2( 1 _ XZ )2 

— z(x — z) 

n _ s 2U/2| 1 _ z 2|l/2n _ x ^) 2 . 



(51) 



Kcc(Zi,Z2) l _ 

" - zg) ~ fg) " 

_1 |l-Z 2 |V 2 |l-z 2 | 1 /2(l_ Zl ^ 2 )2 |l- Z 2 |l/2|l_ z 2|l/2( 1 _ Zl ^ 2 )2 

7T (^1 ~ Z 2 ) i(zi ~ Z 2 ) ' 

|1 - Z 2 |V 2 |l - Z 2 \ l / 2 (1 _ fl ^)2 |1 _ z 2|l/2|! _ 4|1/2(1 _ ^ lZ2 )2 

(52) 

Recalling (f3T|) . we check from these explicit forms that the results ([Sj) , (fTU|) 
for the one and two point real correlations and the result (|16p for the one point 
complex correlation are reclaimed. These results were obtained from the forms ([9]) 
and ()13p . different from the corresponding Pfafnan formula ()30p . An outstanding 
task then is to give a direct derivation of the equivalence between the forms. In 
the case of complex limiting Kac polynomials, the analogous problem was solved 
by Peres and Virag |T7] . 
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